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ABSTRACT 


The importance of supplying undisturbed electricity keep increasing due to 
modernization and lifestyle. Any disturbance in the power system may lead 
to discontinuation and degradation in the power quality. Therefore, detecting 
fault, fault type and fault location is a major issue in power transmission 
system in order to ensure reliable power delivery system. This paper will 


compare two prominent methods to estimate the fault location of double 

circuit transmission line. Those methods are Discrete Wavelet Transform 
Keywords: algorithm and Fast Fourier Transform algorithm. Simulations has been 
carried out in MATLAB/Simulink and a variety of fault has been imposed in 
order to analyse the capability and accuracy of the fault location detection 
algorithm. Results obtained portrayed that both algorithms provide good 
performance in estimating the fault location. However, the maximum 
percentage error produced by the Discrete Wavelet Transform is only 0.25%, 
0.6% lower than maximum error produces by Fast Fourier Transform 
algorithm. As a conclusion, Discrete Wavelet Transform possesses better 
capability to estimate fault location as compared to Fast Fourier 
Transform algorithm. 
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1, INTRODUCTION 

Transmission line is one of the main components in electrical power system which works as 
medium for transmitting electrical power generated to the consumers. In order to ensure a continuous 
electricity supply at the consumer side, maintenance and protection should be in high concern. However, 
most of the time, fault occur without a warning and can be at any place due to weather condition, failure of 
the equipment, error by human mistakes and others. The past two decades has viewed the rapid development 
in techniques for detection, classification and location of faults in power systems. The advancement in 
communication system, signal processing techniques, global positioning system (GPS), artificial intelligence 
system and machine learning have enabled more thorough research to be carried out in fault protection 
techniques. In the past, the locations of fault are hardly to determine due to limitation of technology and 
method in detecting and locating it in the real time monitoring. A key point for a fault location was an 
accuracy and reliability for the algorithm proposed. Fault location methods can generally be classified into 
two fundamental categories: techniques based on power-frequency components, and the other utilizing the 
higher-frequency components of the transient fault signals [1]. 


Journal homepage: http://iaescore.com/journals/index.php/ijeecs 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 1357 


Theoretically, the frequency features of voltage and current will change dramatically whenever a 
fault occurs. If these features able to be identified and interpreted correctly, it may help protect the affected 
power system network to a great extent. It 1s also essential to classify the faults and identify which phases are 
involved. Various methods have been proposed to detect, classify, or locate faults in power transmission lines 
in recent years. One of it is Fourier Transform (FT). FT is a widely used mathematical tool which mainly 
used to analyse the frequency spectrum of a signal. Discrete FT (DFT) is applied when the signal is in 
discrete form. Authors in [2] used half cycle DFT to filter out the harmonic component in the signals and 
estimate its phasor element. [3-5] also adopted half cycle DFT to calculate fundamental and harmonic 
phasors of the signal to classify the fault-type in the studied system. WT is one of the most popular feature 
extraction methods for various fault diagnosis systems. In practice, most of the studies use Discrete Wavelet 
Transform to decompose the original current and voltage signals into its time and frequency spectrum. The 
WT decomposes signals over dilated and translated wavelet, which these compositions known as multi- 
resolution analysis (MRA). In this approach, the high-frequency components of the signal will have a larger 
time resolution while its low-frequency components will have larger frequency resolution [6]. Wavelets have 
shown to perform better with non-periodic signals that contain short duration impulse components as is 
typical in power system transients [7]. 

A study on WT locating a fault has been done by many researches due to the strongest tool to 
analyse the fault because of its perfect time frequency localization ability [8-12]. WT has an advantages in 
detecting fault in transmission line where the WT provided a varies size analysis window in proportion to the 
frequency analysis hence offer a better compromise in terms of localization [13]. A research on FFT to 
locating fault by [14] state that the frequency of the first harmonic is utilized for determination of the fault 
location using travelling wave theory and the results show that the accuracy of the method is not affected by 
the fault resistance and phase angle of fault instant. 

In [16], a classification algorithm of different faults on transmission line using fuzzy logic has been 
presented. FFT was applied in the study to determine the fundamental component of the current signals. 
Reference [18] has compared the performance of FT method with wavelet in detection and classification of 
faults on transmission line. But then, it is reported that wavelet give a significant performance for fault with 
more than single phase involved. Moreover, most fault studied in the paper are single phase to ground faults. 
While in [17], a neural network based fault identification and classification scheme is proposed. Reference 
[19] also discussed the application of neural network in double end fed transmission line for faulty phase 
selection and fault distance location. Authors in [20] has proposed to classify only ground faults on double 
circuit transmission line using artificial neural network. Application of wavelet in simultaneous fault in 
transmission system identification is discussed in [15]. However, authors have suggested to improve the 
proposed algorithm in order to achieve the desired accuracy. 

With the evolvement of wavelet, some researchers proposed to use new type of wavelet to obtained 
useful information from the voltage and current signals for fault identification and classification. In [21] and 
[22], the entropy (energy) of the signal was extracted using Discrete Wavelet Transform and fed to artificial 
neural network (ANN) for fault classification. In [23], the current and voltage signals at relay location is used 
and transformed into series of details and coefficient using Discrete Wavelet Transform. These details and 
coefficient are then fed to ANN for fault detection and classification 

The work presented in this paper deals with fault location detection using Discrete Wavelet 
Transform for three types of faults in double-circuit transmission lines. A 220 kV line double end fed 200 km 
in length double circuit transmission line has been chosen as the studied system. The work reports the result 
of an extensive ‘offline’ studies using Matlab and its associated toolboxes: Simulink, SymPowerSystems and 
Wavelet toolbox. All the three types of faults (single-line to ground fault, line-to-line fault and dounle-line to 
ground fault) has been assigned on both lines in double-circuit transmission line at the same time. The 
aptness and the effectiveness of the selected technique is analysed and compared with the percentage error of 
the fault location using Fast Fourier Transform. 


2. DISCRETE WAVELET TRANSFORM FAULT DETECTION AND FAULT LOCATION 
ALGORITHM 
Figure | below provides the algorithm of implementing Discrete Wavelet Transform which has been 
proposed in [24], in order to detect the fault occurrence on double circuit transmission line. In detecting the 
fault on transmission line, three phase current signal will be recorded, and any inception of fault will change 
the amplitude of the three-phase current signal will be extracted to the wavelet decomposition. Figure 2 
depict the wavelet decomposition process. 
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Figure 1. An algorithm in detecting fault using discrete wavelet transform 
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Figure 2. Wavelet decomposition process 


Wavelet decomposition gains an input from the current measurement in bus terminal. The signal 


then, started to decompose into three levels of decomposition by passing through a low pass filter. The result 
from the decomposition using discrete wavelet decomposition filter will determine the fault inception by 
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obtaining the high frequency detail. If the resulted absolute value from the first difference of high frequency 
detail coefficient is greater than the threshold value, it denotes that fault is established within the transmission 
lines. Once fault detected, the location of the fault will be determined further using: 


FLRTXTL 
FLptT X FLst 


Fault Location: 


(1) 


Where FLST is a single faulted phase from the source terminal while FLRT depict local fault 
location information from the remote terminal and the length of transmission is represented as TL. 


3. FAST FOURIER TRANSFORM FAULT LOCATION ALGORITHM 

The proposed method [25] possess some advantages where fault resistance and in feed from both 
substations do not affect fault location estimation. Also, the different source impedance values will not 
influence the fault location estimation and the data is not required. During fault, the harmonic component was 
added in the voltage and current signal. Thus, an algorithm based on the harmonic component has been 
developed to locate the fault. The faulted voltage and current signal has to be extracted from the transmission 
line model before filtering out the harmonic component from the signal. Figure below shows a voltage and 
current signal extraction. Later, The FFT block functioned to convert the sampled voltage and current 
waveforms into magnitude and phase angle at every samples. Whilst, the low pass filter block functioned as 
filtering a harmonic component. Figure 3 provides the block diagram of fault location using Fast Fourier 
transform. 
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Figure 3. Fault location using Fast Fourier transform 


The output from the FFT block was being converted into a magnitude and phase before the 
estimation of fault location. The fault location later can be estimated using equation given in [25], which can 
be summarized as: 


Va-Vb+ CAI +Z2 +Zo) 


ea (2) 
140 (Z1+Z2+Z) 


Where; 

m = line length (in per unit) 

Z1 = Positive sequence impedance 

Z2= Negative sequence impedance 

Zo = Zero sequence impedance 

Va= Phase to ground voltage seen from substation A 
Vs= Phase to ground voltage seen from substation B 
I, = Phase current from substation A 

Ip = Phase current from substation B 
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4. DOUBLE LINE TRANSMISSION MODEL 

In order to investigate the applicability of the proposed method, a model of double line transmission 
system as shown in Figure 4 is considered for the purpose. The system testing model has been developed 
using MATLAB/Simulink Software. The model is composed by 220 kV line with 200 km in length where 
single source and the short circuit capacity has been set to 1.25 GVA. Xs/Rs ratio of the source is set to 10 
[26]. The parameters of the double line model used in this circuit are as listed in Table 1. 
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Figure 4. Effects of selecting different switching under dynamic condition 


Table 1. Parameter used in the Double Line Transmission Model 


Parameter Value 
Positive sequence resistance, R, 0.01809 Q /km 
Zero sequence resistance, Ro 0.2188 Q/km 
Zero sequence mutual resistance, Rom 0.20052 Q/km 
Positive sequence inductance, L 0.00092974 H/ km 
Zero sequence inductance, Lo 0.0032829 H/ km 
Zero sequence mutual inductance, Lom 0.0020802 H/ km 
Positive sequence capacitance, C, 1.2571e-008 F/ km 
Zero sequence capacitance, Co 7.8555e-009 F/ km 
Zero sequence mutual capacitance, Com -2.0444e-009 F/ km 


5. RESULTS AND ANALYSIS 
5.1. Fault Detection and Fault Location Estimation using Discrete Wavelet Transform 

Various types of faults including single line to ground fault, line to line fault and double line to 
ground fault, has been initiated into the system randomly along the transmission lines. Figures 5, 6 and 7 
illustrates the wavelet decomposition of the current signal during faults. For all types of faults, the fault 
existence has been confirmed after the 3rd level decomposition process via dl value in the decomposition 
equation. Table 2 tabulate the results for faults detection status from using Discrete Wavelet 
Transform algorithm. 


Table 2. Discrete Wavelet Transform Fault Detection in Transmission Line 


Location fault (km) Single line to ground fault Line to line fault Double line to ground fault 
40 Fault detected Fault detected Fault detected 
80 Fault detected Fault detected Fault detected 
120 Fault detected Fault detected Fault detected 
160 Fault detected Fault detected Fault detected 
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Figure 5. A wavelet decomposition of current signal during single line to ground fault 
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Figure 6. A wavelet decomposition of current signal during line to line fault 
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Figure 7. A wavelet decomposition of current signal during double line to ground fault 


From the results, it can be concluded that Discrete Wavelet Transform has the ability in detecting 
fault within transmission line in any cases of fault and without concerning the location of the fault itself. 
Further, Table 3 list down the estimated faults location under various positions along the transmission line. 
Also, only single line to ground fault 1s being considered during the simulation. The highest error recorded 
is 0.5 km 
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Table 3. Estimated Fault Location in Transmission Line using Discrete Wavelet Transform 


5.2. Fault Location Estimation using Fast Fourier Transform 


Actual fault location (km) 


Estimated fault location (km) 


19.95 
39.90 
59.85 
79.80 
99.75 
119.70 
139.65 
159.60 
L955 
199.50 


Error (km) 


0.05 
0.10 
0.15 
0.20 
0.25 
0.30 
0.35 
0.40 
0.45 
0.50 


Whilst taking red phase as reference and implying single line to ground fault towards the double line 
transmission system, Table 4 tabulated the resulted estimated faults location under various positions along 
the transmission line by using Fast Fourier Transform algorithm. Error obtains ranging from lowest at 0.1 km 
to the highest at 1.7km. 


Table 4. Estimated Fault Location in Transmission Line using Fast Fourier Transform Algorithm 


5.3. Fault Location Estimation Comparison 


Actual fault location (km) 


Estimated fault location (km) 


21.2 
40.9 
60.7 
80.4 
100.2 
119.9 
139.5 
158.8 
178.5 
198.3 


Error (km) 


1.2 
0.9 
0.7 
0.4 
0.2 
0.1 
0.5 
1.2 
1s) 
1.7 


Table 5 provide the comparison of estimated fault detection location for single line to ground fault 
between Fast Fourier Transform algorithm and Discrete Wavelet Transform algorithm. Figure 8 and Figure 9 
provide the comparison between both methods in terms of its estimated fault location and percentage errors 
of estimated fault location. In estimating the fault location, the maximum magnitude of the percentage error 
recorded is 0.85%, which is obtained from Fast Fourier Transform. But the magnitude of error produced by 
implementing Discrete Wavelet Transform algorithm is mainly lower compared to the one produces from 
Fast Fourier Transform algorithm. The highest percentage error obtained by Discrete Wavelet Transform is 
only 0.25% where actual fault location is set to 200km. 


Table 5. Estimated Fault Location in Transmission Line using Fast Fourier Transform Algorithm and 
Discrete Wavelet Transform Algorithm 
Discrete Wavelet Transform 


Actual fault location 
(km) 


20 
40 
60 
80 
100 
120 
140 
160 
180 
200 
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Estimated fault location 
(km) 
19.95 
39.90 
59.85 
79.80 
99.75 
119.70 
139.65 
159.60 
179.55 

199.50 


Error (%) 


-0.025 
-0.050 
-0.075 
-0.100 
-0.125 
-0.150 
-0.175 
-0.200 
-0.225 
-0.250 


Fast Fourier Transform 


Estimated fault location 


(km) 
65 
40.9 
60.7 
80.4 
100.2 
119.9 
139.5 
158.8 
178.5 
198.3 


Error (%) 


0.60 
0.45 
0.35 
0.20 
0.10 
-0.05 
-0.25 
-0.60 
-0.75 
-0.85 
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Figure 8. Estimated fault location in transmission line using Fast Fourier Transform algorithm and Discrete 
Wavelet Transform algorithm 
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Figure 9. Error of estimated fault location in transmission line using Fast Fourier Transform algorithm and 
Discrete Wavelet Transform algorithm 


6. CONCLUSION 

The proposed method utilizes Discrete Wavelet Transform spectral energy values for fault location 
identification. Performance comparison of the studied method has been performed with Fast Fourier 
Transform. Both of the algorithms portray good performance of identifying fault location. However, results 
obtained indicate that the Discrete Wavelet Transform demonstrates higher accuracy capability as compared 
to Fast Fourier transform. In addition, Discrete Wavelet Transform algorithm also be able to detect the fault 
occurrences apart from estimating the faults location. 
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